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Abstract 

We analyze CP asymmetries of B — > <j)Ks and .B — ► n'Ks in a supersymmetric 
grand unified theory in which only the third generation sfermions contained in 
10(Q,U C ,E C ) of SU(5) can have a different mass from the others. One of the 
advantages of this nonuniversal mass model is that the first two generation sfermion 
masses can be large whereas both (left and right handed) stops are light so as to 
stabilize the weak scale. Therefore, we studied a minimal supersymmetric standard 
model parameter region in which a fine tuning in Higgs sector is relaxed owing 
to light masses of stops, gluino and higgsinos. In such a parameter region, the 
chargino contribution is as important as the gluino one. We show that the CP 
asymmetries of B — > (pKs and B — > n'Ks can deviate from their standard model 
predicted values by 0(0.1) because of constructive interference between gluino and 
chargino contributions. 



1 Introduction 

Super symmetry (SUSY) is one of the most promising candidates for physics beyond the 
Standard Model (SM). The minimal supersymmetric standard model (MSSM) not only 
provides a solution to the gauge hierarchy problem in the SM but also has some attractive 
features, for example, the gauge coupling unification, the radiative electroweak symmetry 
breaking, and giving dark matter candidates as the lightest superparticle. 

However, the MSSM has serious problems in SUSY breaking sector. If we introduce 
generic SUSY breaking terms, they induce very large amplitudes for Flavor Changing 
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Neutral Current (FCNC) and CP violating processes to satisfy the various experimental 
constraints [H, El [3]. To avoid the problems the universal soft sfermion masses are often 
assumed at some scale [U [5]. 

The assumption of the universal soft masses is not necessarily required. For example, 
since the FCNC constraints for the first two generation field are much severer than for the 
third generation fields, the universal soft masses sometimes have been imposed only for the 
first two generation sfermion, which are realized if we impose non-Abelian horizontal sym- 
metry, for example, U(2), under which the first two generation fields are doublets and the 
third generation fields are singlets [7j. However, if the diagonalizing matrix of the fermion 
has Maki-Nakagawa-Sakata [S] (MNS)-like large mixings, such sfermion mass spectrum 
leads to very large FCNC to be consistent with the experimental bound. Therefore, such 
nonuniversality should be introduced only for the sfermion whose fermionic superpartners 
have the diagnalizing matrix with the Cabibbo-Kobayashi-Maskawa [9] (CKM)-like small 
mixings. 

In SU(5) SUSY Grand Unified Theory (GUT), it is a reasonable assumption that 
the diagonalizing matrices of 10 have the CKM-like small mixings, while those of 5 have 
the MNS-like large mixing, because the 10 includes the doublet quarks, and 5 includes 
the doublet leptons. Under this assumption, only the sfermions included in the third 
generation of the 10 (IO3), can have different masses from the others without conflicting 
with experimental constraints from various FCNC processes. 

We summarize several characteristic features of the nonuniversal sfermion masses. 
First of all, the rotation matrices for 10, which make fermions to mass eigenstates are 
expected to have the CKM-like small mixings, since 10 involves quark doublet. Therefore, 
large off-diagonal entries for the 10 sfermion mass matrices do not arise after this rotation, 
even if the initial soft masses are not universal. Moreover, FCNC constraints among 1-3 
or 2-3 generations are not so severe compared with those of 1-2 generations. Therefore, 
we can expect that nonuniversality for the third generation does not conflict with the 
FCNC constraints. Moreover, it is important that we can expect larger FCNC for the 
third generation fields than in the usual universal sfermion mass case as discussed in 
Ref . [T£] . Second, the naturalness of the Higgs mass in MSSM requires that the gaugino 
masses, the higgsino mass, and the stop masses must be around the weak scale. In the 
nonuniversality case, we can take the larger masses for all the sfermions except for the 
sfermions of the IO3 without conflicting with the naturalness arguments, because both 
(left and right handed) stops are included in the IO3. In principle, we can take such 
mass larger than 1 TeV, which can relax the various constraints from FCNC and electric 
dipole moments (EDM), etc. Again, the FCNC related with the third generation fields 
can become relatively large, which may be detectable in future experiments. 

Finally, such a sfermion mass spectrum can be realized in the anomalous U(l) SUSY- 
breaking models[10j, flavor-mediated SUSY breaking modelspT], radiatively driven mod- 
els with specific boundary conditions [T2], extra dimension models, [TS] and models with 
discrete symmetry, [TJ] and is naturally derived from E 6 SUSY GUT model [To"| ITo] with 
SU{2) or £77(3) horizontal symmetry [T7] . 

It is expected that FCNC processes among 2-3 generations, particularly, the bi — > sl 
or tr — > fift transition rate become relatively large in the nonuniversal sfermion mass 
model. In the previous paper [IB] , we studied r — > ^7 and /i — > e-y processes and 
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showed that these lepton flavor violating processes can have very large branching ratios 
and that these processes may be discovered in near future experiments. However, in 
the quark sector, this model gives similar predictions for the various FCNC processes to 
those in the model with the universal sfermion masses because the off-diagonal elements 
of the squark mass matrix obtained from the nonuniversal sfermion masses are of the 
same order as those obtained by the renormalization group equations (RGEs) from the 
universal sfermion masses. The difference between these models appears in the CP 
violating processes because, in the nonuniversal sfermion mass model, there are new CP 
phases except for the Kobayashi-Maskawa (KM) phase. Therefore, it is interesting to 
study the CP violating processes in the nonuniversal sfermion mass model to distinguish 
these models. 

In this paper, we focus on the 2-3 transition in the quark sector and analyze the CP 
asymmetries of B — > <f)Ks and B — > rj'Ks [19]. These observables attracted attention 
a few years ago since the Belle Collaboration [2U] reported the large deviation from a 
SM expectation sin 20! (sin 2/3) = 0.681 ± 0.025 [21]. Their current experimental world 
averages are S^k = 0.39 ±0.17, S^k = 0.61 ±0.07 [2H 22j and they are almost consistent 
with the SM expectation, but it is still important to measure these observables accurately 
in the search for new physics beyond the SM. Particularly, the CP asymmetries of B — > 
4>Ks and B — > rj'Ks are interesting for the model with the nonuniversal masses because of 
the following reasons, i) Since SUSY contribution to CP violation in — > transition 
can be large in this model while the SM has little CP phase in the b — > s processes, 
large deviations of CP asymmetries from the SM are expected, ii) Both experimental 
central values of B — > <pKs and B —>■ rj'Ks negatively deviate from the SM expectation. 
It is known that such deviations favor the models with only LL mixing (not RR mixing) 
in down-type squark sector [23]. This gives a strong motivation to analyze the CP 
asymmetries of B — > <fiKg and B — > rj'Ks hi this nonuniversal sfermion mass model, 
because the 10 contain left-handed down quarks and not right-handed down quarks, 
iii) As in the processes of r — > ^7 and ri — > ey, the deviations of CP asymmetries 
of B — > 4>Ks and B — > rj'Ks from the SM expectation strongly depend on the sfermion 
masses of the IO3. Therefore, if the sfermion masses of the IO3 are measured in the Large 
Hadron Collider (LHC) experiment, the predictions of the branching ratios of the LFV 
processes and the CP asymmetries can be more definite. Even if the sfermion masses 
of the IO3 are not measured by LHC, we can check the consistency of this model by 
comparing with the LFV processes and the CP asymmetries. 

This paper is organized as follows. In the next section, we will specify the model and 
discuss the flavor and CP violating couplings in this model. After giving definitions and 
notations to analyze CP asymmetries in §3, the results of numerical study are given in 
§4. In § 5 we will discuss some of the constraints of this model. We summarize this paper 
in §6. 
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2 CP violation in the nonuniversal sfermion mass 
model 



In this section, we discuss flavor and CP violating couplings related with CP asymmetries 
of b — > s processes in the nonuniversal mass model. First, we assume that the SM 
gauge group SU(3)c X SU(2)l x U(1)y is unified into SU(5) and the mediation of SUSY 
breaking respects the SU (5) gauge symmetry. We start from the following Lagrangian 
at GUT scale, 

C = u R ^q j H u + d Ri rf j ' e q j H d + --- 

+ u* m (ml )iiu m + q*(m 2 w )uqi + d Ri (ml)ud m + • • • , (1) 

where q iy u c Ri , d c Ri (i = 1,2,3) are quark doublets, right-handed up quarks, and right- 
handed down quarks, respectively, and the fields with tildes denote sfermions. As men- 
tioned in the Introduction, we assume the following forms for the sfermion mass matrices. 

(m% \ (ml \ 

m 2 = ml , m\ = m\ , (2) 

V m loJ \ m lJ 

By redefining the superfields, we can diagonalize the up type Yukawa matrix. 

U R ->VuU R , D R ^V D D R , Q^V Q Q. (3) 

C = u R% {Y-) llULl Hl + d R% {% d V^^ 

+ u*BA m l R o)ijURj + qt( m2 qo)m + d Ri (ml)iid m + ■■■ , (4) 

where Y f and V CKM are the real and diagonal Yukawa matrices and CKM matrix at 
GUT scale, respectively. Vq is a unitary matrix of quark doublets and it is roughly 
expected as 

(1 A A3\ 

V Q ~ A 1 A 2 , (5) 
\A 3 A 2 1 / 

where A is Cabibbo angle (A = 0.22) and we omit the 0(1) coefficients and phases. Then, 
m? is estimated as 

'ml 



m 1o = V Q m W V Q ~ m l \+ m 



m FC0 — Vq 






|v« 


ml) 




A 5 


A 3 N 


A 4 


A 2 


A 2 


j 



m 2 - m 2 0/ 



mi, -mi) A 5 A 4 A 2 . (6) 



1 Note that the following arguments can be applied to unified models with larger unified gauge groups, 
for instance, 5*0(10) or Eq. 
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Here, we use m\ = m^l + diag(0, 0, m| — m 2 ,) and the unitarity properties of Vq. Since 
we have already used the freedom of phase rotation of quark fields to reduce the number of 
phases in Vq KM into one, the coefficients in Vq have 0(1) phases generically. Therefore, 
the off-diagonal entries of ra~ and m 2 FC0 can have new 0(1) phases. 

To calculate CP asymmetries of B — > <p(rf)Ks, we need the low energy sfermion mass 
matrices m~ and m 2 ~ , which are obtained by renormalization group equation (RGE). 
However, off-diagonal entries of sfermion mass matrices do not change significantly in 
RGEs between GUT and weak scale, for example, 



\Y U 1 2 1 

(m|)32(/Wfc) ~ exp[ ™ \og(fil eak / A 2 GUT )] x (m|) 32 (A GC7T ) 

~ 0.8 x (m 2 ) 32 {A GUT ) . (7) 

In our analysis, we do not fix the off-diagonal entries at GUT scale (We can estimate only 
their order.), so we absorb the high and low energy differences of the (mj) 3 2 in ambiguity 
of (771~)32(Ag(7t). For other relevant off-diagonal entries, the changes between GUT and 
weak scale are rather small owing to small Yukawa couplings. 

Therefore, we calculate the diagonal components of these mass matrices by RGEs. 
Finally, we go to the super CKM basis by rotating the superfields a.s D L ^ VckmDl at 
low energy. Then, low energy sfermion mass matrices are obtained as 

'm 2 \ 

ml L ~ ( m 2 J + m 2 FC , (8) 



m 3, 



2 

m 



m 2 r , T/t I ™2 



v ckm m2 Vckm + V CKM m 2 FC V C KM, (9) 



m 2 ~ does not have off-diagonal entries. Note that since a component {Vckm)z2 does not 
have a large CP phase, the first term in Eq.Q does not contribute to the parameter 
Im[(m| ) 32 ]. Thus, the following approximate relation is obtained. 

Im[(my 32 ] ~ Im[(m?J 32 ] ~ (m 2 30 - m 2 )A 2 sin0 SUSY , (10) 

where S usy is a phase of (mpc) 32- This relation means that if the gluino contribution is 
maximized by the SUSY phase, then the chargino contribution is also maximized. Thus, 
strong interference between gluino and chargino contributions is expected in this model. 



3 CP asymmetry of B — > (f>Ks and B — > rj'Ks 

In this section, we give a review of the well known method of effective field theory for- 
malism to calculate the CP asymmetries of B — > (f)(j]')Ks [2^1 125]. Their time dependent 
CP asymmetries are defined as 

T(B(t) -> tfjQKs) - T(B(t) -> tfjQKs) 
a ^')KsW - r(B(t)^mKs)+r(B(t)^^rj')K s ) 

= A%, )Ks cos AM Bd t + StMKs sin AM Bd t, (11) 
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where A^L,\ Ks and S^^Ks represent the direct and mixing CP asymmetries, respectively, 
and their expressions are given as 

dir |p(0(r/)if s )| 2 -l 2Im[e~ 2 ^p(0(r/)iG)] 

where 0i(or 0) is the standard angle of the unitarity triangle, and parameter p(<fi(r)')Ks) 
is defined as 

p(0(^)^) = A(W) ^ r (13) 

Here, A((f>(r]')Ks) and A(cf)(rj')Ks) are decay amplitudes of P — > (p(rj')Ks and P — > 
(j>(r/)Ks, respectively, which can be calculated using the effective Hamiltonian of AB = 1 
transition at the low energy fi b ~ m b 



A(4>( V ')Ks) = (mK s \H^ =1 \B), A(4>( V ')K S ) = {4,{ff)K a \{H^=^\B). (14) 

tAB- 
l eff 



l eff I-"/ > -n-KVVl )^S) — WVU-^SWeff 

Effective Hamiltonian PfV? =1 is expressed using the operator product expansion(OPE) 



as 

6,7 7 ,8 9 

^e A /f =1 = E E S t^^ 1 ^ + - + W - 4 C - C}, 

p=u,c i=3 

(15) 

where X q = V q bV* s , Ci = Cj(/ib) are the low energy Wilson coefficients. The low energy 
renormalized operators = Qi(//&) are expressed as 

QS P) = (p a ^2P L b p )(s^2P LPa ), = (p lfl 2P L b)(sY2PLP), 

Q 3 = (s lfi 2P L b) ^(g 7 ^2P L g), Q 4 = {s al>J 2P L bp) ^(^2P ife ), 

Q 5 = (s 7m 2P l 6) ^(g 7 ^2P fl g), Q 6 = (s a7fl 2P L bf,) E^T^fe), 

Gr 7 = ^m b sa^2P R F^b, Q Sg = ^m b s a a^2P R G^T*b p , (16) 

where a and /3 are color indices, is SU(3)c generator, a^ u = | [7^,7^], Pl,.r = (1 T 
7s)/2 represent the projection operators, and Qj is obtained from Qi by exchanging L <->• 
P. Here, we ignored the electroweak penguin operators Q7-10 an d the contributions to 
dipole operators which are proportional to the strange quark mass. The matrix elements 
(B(B)\Qi\(j)(ri')Ks) are provided in Appendix A. The low energy Wilson coefficients Cj 
can be obtained from the high energy coefficients Ci(nw) (pw — m w) by solving the 
renormalization group equations for QCD in the SM. The solution is expressed as [24] 



Ci(Pb) = E^'(^6' Pw)Cj(nw)- (17) 



Here, Uij is the evolution matrix and we included the leading order in QCD. Their ex- 
pressions are provided in Appendix B. For simplicity, the matching scale is chosen as the 
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raw, but we calculate the Wilson coefficients in full theory, which includes the super- 
particles. We construct each coefficient using the SM contributions Cf M (iiw)i Charged 
Higgs contributions C^fiw), Gluino contributions Cf {nw)i an d Chargino contributions 
Cf(iA-w) as follows: 





(jSM 






c 2 = 


r>SM 




c 3 = 


cl M + ci 






C A = 


^-(SM j_ n g 




c 5 = 


c! M + cf, 






c 6 = 


°6 ' 6 ' 




C77 - 


r<SM , 
- c- ?7 -t- o 77 


+ cf 7 






^8g + °8g 


+ C 8g + Cg fi 



Here, we ignored the contributions of C| / (/j,w)-Cq \hw) an d ^ipwY^ipw)- H Since 
there are no flavor and CP violation in right-handed down- type squark sector in the 
nonuniversal mass model as mentioned in the previous section, the coefficient C 9 does 
not emerge. The C H and C x can also be neglected because the CV 7i $ g is suppressed by the 
factor m s /rrib compared with Cr^sg- The difference in the final state parity of B — > (pKs 
and B — > rj'Ks leads to the following relation. 

(0#°|Q<|fl°) = (^UlB ), (rj'K'mB ) = -WK'UIB ). (19) 

The 1-loop order expression for the Wilson coefficients in Eq. fllSI) are provided in Ap- 
pendix B. 

Let us discuss the deviations of Sqk an d S^k from the SM expectation. First, we 
extract the strong phase from the amplitudes as A = e tSs A CP . Since the dominant 
contribution to the A((f>(rj') Ks) , which is provided from the SM is almost real, we can 
expand p{(j){i)K s ) and p{<f>{rf)K 3 ) with R^ )Ks = lmA CF \<j>{rf)K s ) /ReA CF \<j>{rf)K s ) 
as follows. 

ReA CP - ilmA CP , . . 

'- w + anAo * 1 -**-**- (20) 

Then using Eq. (fT2l . we can derive 

S<f>(n')K s = sin20x + AS , 0( r? /)^ s , (21) 
ASVw)^s = 2[R^ v > )Ks cos2(j) 1 - Rl w)Ks sin2<j)i], (22) 

where sin20x = 0.68, cos 20! = —0.73, and AS^^Ks represent a deviation from the 
SM expectation. The first term provides the leading contribution to AS^^Ks, and 
second term always provides the negative contribution. Therefore, the maximum value 
of \AS^ v i)k s \ is given as positive R<f>(ri')K s - 



4 Numerical analysis 

In this section, we show the numerical results of the deviation of CP asymmetries of 
B — > 4>Ks and B — > rj'Ks from the SM expectation. We are interested in the parameter 
region in which a tuning in the MSSM Higgs sector is relaxed owing to the light masses 
of stops, gluino and higgsinos, because one of the advantages of the nonuniversal mass 

2 We checked numerically that their contributions cannot be large in our parameter region of interest. 
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model is that the light stops can be realized while the first two generation sfermions are 
heavy. Therefore, we basically use the following parameter set through our analysis. 

At GUT scale : 

m 30 = lOOGeV, m = 1000 GeV, m 1/2 = 200 GeV, A u = A d = GeV; 
at weak scale : 

fj, = 250 GeV, tan/3 = 10, 9 SUSY = -tt/2, (23) 

where mu^ is the universal gaugino mass and /i is the higgsino mass. We took them as 
relatively small to relax the fine tuning in the MSSM Higgs sector. 

At first glance, we may think that the above parameter region conflicts with the LEP 
constraint for the higgs boson mass. However, the well-known LEP constraint on the SM 
Higgs boson mass, 114.4 GeV cannot be applied to the "MSSM" lightest higgs boson if 
the MSSM lightest higgs does not have relevant coupling to the Z boson [261 E3 EE]. In 
this situation, the large quantum correction to the SM-like Higgs boson from the large 
stop and gluino masses is not required. We checked that the above parameter region 
is numerically consistent with the LEP constraint in the literature In the above 

parameter region, since charginos and stops are relatively light, it is expected that not 
only gluino diagram but also chargino diagrams give relatively large contribution to the 
CP asymmetries. 

In such a small ZZh coupling region, it is known that all MSSM Higgs bosons have 
relatively small masses [2"fi"| |2"T] . Then, negative /i and large tan/? are disfavored by the 
constraints from b —> S'j, t — ► H + b, B£ — > r + u T , and B s — > processes. There 

are large contributions to the h — > 57 from the charged Higgs-top loop diagram and the 
chargino-stop loop diagram in our parameter region. The typical orders of magnitudes 
of their amplitudes are of the same order as the SM one and the cancelation between 
them is therefore required. This cancelation is realized only in the positive \x case. At 
the same time, tan/3 cannot be arbitrarily large otherwise the top quark significantly 
decays into charged Higgs boson and b quark (tan/3 < 40 is mandatory [29] ). Moreover, 
the contributions from the light Higgs boson exchange to the BR(B^ — > t + v t ) and the 
BR(B S — > are highly sensitive to large tan/3, because they are proportional to 

(tan/3) 4 and (tan/3) 6 , respectively [3U| I3"T] Thus, the constraints from these processes 
disfavor the region in which tan/3 > 15[32j. From the above reasons, we do not consider 
the negative \i or large tan (3 region in this paper. 

Before presenting the numerical results, we give a rough argument to show that we 
can discuss the magnitude of each contribution to AS separately. To separate each 
contribution, we use the following notation: 

A°r{<M)K 3 ) = A$$ )Ka + Afa )Ka + A^ )Kg + A<* )Ka + A*, )Ka , (24) 



8 



where A^,^ Ks is the SM amplitude and the remainders are 



A Hv')k s = -^\t[0.727C^ w )m V ')Ks\Qs 9 \B°), 



AX 

A 4>{rf) 



V2 
G F 



Ks = -^A i [O.727^(^)](0(r 7 ')^|Q 8s |S ), 



Gi " " 



A MKa = --/§ A * X)E &i(^> Vw)C?(fi w )] W)K s \Qi\B°), 
»=i i=i 

= -%A 4 [O.727Cl q ( / i W )](0(r ? ')^5|Q 89 |5 ), (25) 



where the factor 0.727 can be found in (1561) in Appendix B. Noting that the imaginary 
part can appear only in the SUSY contributions, we expand Rm^Ks up to the leading 
order of Re[A x }/A SM (X = H, X , #36, g8) as follows: 



lm[A] 

A H Re [A*] Re [A 936 ] Re[^ 8 ]^ 

(27) 



ASM y[SM asm a sm 

In our interesting parameter region, T ~ 1 because RefA 936 ]/^ 5 *' 1 , Re[A 98 ]/A SM <C 1, and 
A H / A SM ~ — Re[A x ]/74 5M ~ (9(0.1) numerically. Therefore, we obtain the expression of 
R<t>(r)')K s up to leading order as 



ASM 1 ASM 1 ASM 



(28) 



Thus, each contribution can be considered separately. Speaking in more detail, we have 
two kinds of chargino contributions. One of them comes from the superpartner of the 
diagram for the charged Higgs contribution in which flavor violation is originated by 
Vckm- This almost dominates Re[A x ], and thus, we have the above relation |Re[74 x ]| ~ 
|Re[A H ]| in our parameter region because all the mass scales in the chargino diagram are 
of the same order as the mass scales in the charged Higgs diagram. However, this gives 
little contribution to Im[A x ] because Im^VcifAf^] <C Re[(VcfK-Af)32]- The other chargino 
contribution is from (m| L ) 32 , which is caused by the nonuniversality of sfermion masses. 
This gives large contribution to Im[A x ]. 

We present our numerical results for mo dependence of AS^k s an d AS v >k s in Fig. 1. 
In the numerical calculation, we use AS£^ K = —lA6Im[A*,^ K ]/A^£-. K . To make the 
result more correct, for the total ASm^k in Fig. 1, we use Eq. (fT2l . although the sum 
of each contribution is almost equal to the total AS^^k- In the upper figures in Fig. 1, 
the Cg contribution, the Cf gluino contribution, and the Cf-Cf gluino contribution are 
denoted by red (deep gray), green (middle gray), and cyan (light gray) lines, respectively. 
The SUSY contributions are proportional to the sin. 9 susy- Hence, we take 9susy to 
be — 7r/2 through our analysis so that the SUSY contributions are maximized. In our 
numerical calculation, we scan the order one parameter and phases in the matrix Vq in 
the range (0.8 - 1.2) and (0 - 27r) with fixing the Osusy- 
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Figure 1 shows that if mo is much larger than 77130, the CP asymmetries of B — > (pKs 
and B — ► rj'Ks deviate from the SM expectation. At the point mo = m3o = lOOGeV, 
since CP violating couplings (m| ) 32 and (m| i ) 32 vanish, AS^Ks anc ^ ^^rfKs become the 
values predicted by the SM. S^Ks and S^Ks have a SM contribution that originated from 
the Kobayashi-Maskawa phase. From our approximation, the SM predictions of AS^Ks 
and AS v 'k s are 



AS^KsiSM) = +0.02, A3^ Kg (SM) = +0.006. 



(29) 



These SM contributions are indicated by horizontal lines in the upper figures in Fig. 1. 
Figure 1 shows that AS^Ks an< ^ ^^rfK s have the same sign. This is a consequence of 
absence of the RR mixing in the down-type squark sector. As mentioned in § 2, the 
chargino contribution has the same phase as the gluino one. Numerical calculation reveals 
that this interference is constructive. As we can see from Fig. 1, the Cf gluino contribution 
decreases with increasing m , while the Cf-Cf gluino ones and C§ chargino one do not 
decrease. To understand this feature, we use the mass insertion diagrams. The diagrams 
A, B, and C correspond to Cf gluino contribution, Cf-Cf gluino contribution, and Cf 
chargino contribution, respectively. As we can find from the diagrams B and C, for Cf-Cf 
gluino and Cg chargino contribution, cl and §l become heavy with increasing m , but 
this decoupling effect is cancelled by the enhancement of the flavor changing coupling 
(m| J 32 which is proportional to (mj — m^). On the other hand, since the Cf gluino 

contribution has two propagators, which include heavy sfermion (6r and §l) and one 
flavor changing coupling, this contribution is decoupled in the limit m — > 00. 
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Figure 1: too dependence of AS^Ks an d AS 1 ,,/^ and the SUSY contributions. Blue, red, green, and 
cyan (or the order of deeper gray scale) correspond to AS'^/jXsi C*8 chargino contribution, Cf gluino 
contribution and Cf-Cf gluino contribution, respectively. 



We present the numerical results of /i dependence of AS^k s and AS^Ks * n Fig- 2. 
As we can see from Fig. 2, the Cf gluino contribution decrease as decreasing /i, while 
the Cg chargino contribution slightly increase. The reason is that Cf gluino contribu- 
tion (diagram A) is proportional to /i, while C§ chargino contribution (diagram C) is 
decoupled for large \i parameter, because higgsino states are in the loop. The naturalness 
argument requires a small fi parameter; thus, the gluino contribution cannot be large, 
so the chargino contribution and constructive interference between them is important to 
make the sizable deviation of CP asymmetries from the SM expectation. 

Next, we show m 30 and mi/ 2 dependences of AS^Ks anc ^ ^^rfK s 111 Figs. 3 and 4, 
respectively AS^Ks an< ^ AS v ik s strongly depend on m 30 and m^. Again, smallness of 
these parameters is required by the naturalness argument. Thus, the deviations of the 
CP asymmetries from the SM expectation are consistent with the naturalness argument 
in the MSSM Higgs sector. In Fig. 3, we find that the Cg chargino contribution is more 
enhanced than C| gluino one. 

For the chargino contribution, since both stops (t L and in) are included in 10 3 , if ml 
is around the scale A t mt, then one of the mass eigenstates of the stops becomes much 
lighter than the other squarks. As the result, \AS\ ~ 0.1 can be obtained. Again, the 
chargino contribution is important in this parameter region. 

Both the present experimental central values of the CP asymmetries of B — > <pKs and 
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Figure 4: m 1 / 2 dependence of AS<f,K s and AS v 'k s an< i the SUSY contributions. 



B — > r/i^s negatively deviate from the SM expectation as (9(0.1). These are consistent 
with our numerical results. However, we can not extract the new physics contributions 
to | AS 1 1 from the present experimental data because the experimental errors are still 
large and it denotes that the data are almost consistent with the SM predictions. If 
the errors can be reduced in the near future experiments, the deviations from the SM 
can be confirmed and we may obtain several new constrains to our model from the CP 
asymmetries. 

5 Constraints and discussion 

In this section, we discuss the constraints from the b — > sj process and the electric dipole 
moment (EDM). Since CP violating sources for the B — > <pKs and B — » rj'Ks come from 
the flavor changing coupling (m| )32 and (771-^)32 in the nonuniversal mass model, we 
have to explore also the b — > 57 constraint. The experimental value of the branching 
ratio of b -> 57 process Br{b -> s^) exp = (3.55 ± 0.24±°;°? ± 0.03) x 10~ 4 [21] is now 
almost consistent with the SM prediction Br{b — > sj)sm = (2.98 ± 0.26) x 10~ 4 for 
E 1 > 1.6 GeV [33], (3.15 ±0.23) x 10~ 4 for E 1 > 1.6 GeV [M], (3.57 ± 0.49) x 10" 4 for 
E y > m b /20 [33] • Thus, we plot the region in which the total value of the coefficient 
C\ otal {mb) is consistent with the SM value Cj (m&) at 10% level. Figure 5 shows the 
complex plane of C\ otal (m&) / C^ M (m&) . We scan the 9susy in the range (0 - 2tt) in this 
calculation. The grey arrow represents the contribution that includes the SM, charged 
Higgs, and CKM originated chargino contributions, which cannot have the CP violating 
phase except for a small KM phase. As mentioned in the previous section, because of the 
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Figure 5: b — > s-f constraint. We use the parameter set in (ffij)) and m H ± — 130 GeV. 



cancelation between the charged Higgs and the CKM originated chargino contribution in 
the positive \i case, the sum of these contributions becomes almost the same as the SM 
contribution as pointed out in Ref . [26] . On the other hand, the black arrow represents 
the sum of the gluino and {m\ L )^2 originated chargino contributions that have a SUSY 
phase S usy- As we can see from Fig. 5, Osusy — — vr/2 is consistent with the b —> 
constraint. 

For the EDM constraints, it has been pointed out that Chromo EDMs (CEDMs) and 
EDMs strongly constrain the off-diagonal entries of sfermion mass matrices in Ref. [37] . In 
the parameter region adopted in our paper, most of SUSY contributions to the (C)EDMs 
can be decoupled in the large m limit. However, there is a non decoupling diagram con- 
tributing to the (C)EDMs of up quark because both left and right-handed stops have 
small masses. This diagram makes a generically effective mass insertion parameter as 




_«>__^ K> y Im[(^ fl )n]/m^ 

fe 3 io 3 i R 3io 3 u R (™>l L )i3 A t m t (mjjjsi 



2 2 2 2 

m UL m t L m i R 

{m 2 m - mg)A 3 A t m t {m 2 m - mg)A 3 
m 2 m\m\ m 2 

( m 2 )/ m 3 • (30) 

" L 3 



Although accidental cancellation in Im[(5^ R )n] can satisfy these constraints, this is- 
sue can be solved in explicit models. For example, we can construct such models in a 
framework of spontaneous CP violation in E 6 GUT with horizontal symmetry [36J. 

Finally, we comment on the significance of our results for the specific models. In 



3 Actually, exact cancellation is not required here, because there are other contributions for example, 
gluino and (w| )32 originated chargino contributions. 



14 



the specific models, other tree-level contributions to the off-diagonal part of the sfermion 
mass matrices exist generically. Especially, the correction to the 5bar squark fields can 
give a large effect to the result. In this paper, however, we do not discuss them because 
these contributions are dependent on the explicit models and in many cases they can be 
decoupled in the limit m ^> m 3 . If the condition m ^> m 3 is not satisfied it is possible to 
expect more deviation of the CP asymmetry of B — > (f)(i] f )Ks in such explicit models. 

6 Summary 

We investigated the CP asymmetries of B — > <pKs and B — > rj'Ks in a nonuniversal 
mass model in which the third generation sfermion included in 10 of SU(5) have dif- 
ferent masses (~ m^) from the others (~ m). We chose MSSM parameters m ^> m^, 
mi/2 and /i, which suppress the FCNCs and EDMs for the first two generation fields 
without destabilizing the weak scale. Then the chargino contribution is as important as 
the gluino contribution. The reasons are as follows. For the smaller fi parameter, the 
gluino contribution, which is proportional to /itan/3 becomes smaller and the chargino 
contribution becomes larger. Moreover, since both stops are included in 10, small 
makes one of the stops very light owing to the left-right mixing. Owing to the SU(2)l 
symmetry, the chargino contribution has the same CP phase as the gluino contribution 
and the numerical calculation reveal their interference is constructive. As the result, | AS\ 
becomes comparatively large. 

Other CP violating observables related to the b — > s processes, for example, CP 
violation in b — > S7[38] and time dependent CP asymmetry of B s — > J/^0[39], etc., may 
deviate from the SM predictions. This is our future work. 

We found that the deviations of the CP asymmetries of B — > <pKs and B — > f]'Ks 
can be 0(0.1) in this model. We checked that both deviations AS^k and AS v >k have 
the same sign. This is a consequence of which this model has almost only LL mixings. 
These deviations may be confirmed in future experiments. 
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A Hadronic Matrix elements 

Here, we give the expression of hadronic matrix elements for B — > <pKs and B — > rj'Ks 
in the method of the naive factorization. For B — ■> </>lfs, the matrix elements with the 
operator Qi given in ( fl6l) are given as [25j HQ] 



(^IQiIS ) = 0, (<j>K°\Q 2 \B Q ) = 0, (<PK°\Q 3 \B°) = 

(0K°|Q 4 | J B°) = ^X, (^IQsl^ ) = X, (^igel^ ) = (31) 



and 



(^°|q 77 |b°) = o, 

<0X°|Q 89 |5°> = -^_^-[(Q 4 > + (Q 6 > - l((g 3 ) + <Q 5 »], (32) 
where X is defined as 

X = 2F^ K {ml)Um^{p K ■ H ). (33) 

Here, represents the meson mass, is decay constant of the meson, F+~^ K is 
the transition form factor which is estimated at the m$ scale, px is momentum of the K 
meson, and represents the polarization vector of the (ft meson. We used = 1.02 GeV, 
/ = 0.233 GeV, F^ K = 0.35, and 



where mg and rrix are S meson mass and K meson mass, respectively, (q 2 ) in (1321) is the 
average of momentum carried by virtual gluon in Q 8g . Kinematical consideration leads 
to the physical range ml/ 4 < (q 2 ) < m 2 /2. We use (q 2 ) = ml/4 numerically, in this case 
SUSY contributions to CP asymmetries are maximized. Hadronic matrix elements for 
the operator Qi are given by 

(<f>K°\Qi\B ) = (tK^QilB ). (35) 

This relation is derived from the parity invariance of the strong interaction and the fact 
that the initial and final states have the same parity. 

Hadronic matrix elements for B — > rj'Kg are given by [231 HO] 

(iK°\Q x \B ) = X 2 , ( V 'K°\Q 2 \B ) = ±X 2 , 

(r,'K°\Q 3 \B°) = \x x + 2X 2 + ^X 3 , ( V ' K°\Q 4 \B°) =X l + jjx 2 + jjx 3 , 
( V 'K°\Q 5 \B ) = ^X x - 2X 2 - (1 - ^)X 3 , 

(TfK°\Q t \B°) = R X X X - \x 2 - (i - R 2 )X 3 , (36) 
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and 



(r/K u |Q 77 |73 u > = 0, 

(v'k \Qs 9 \b°) = -g-^=[(Q 4 > + (Q 6 > - i«Q 3 > + <Q 6 »], 



(37) 



where 



Xi = -(ro£-mJ)Jf-~( 



m 



2 A^r,' „ 



K 



V2- 



X 2 = -(m%-m$)F*-"(rnh)^fK, 



Xz = -(m z B -m z v ,)Fl 



B^K, 



m: 



Ri 



2m 2 K 



(m b - m s )(m s + m d ) 



2m 2 K — ml 



(m b - m s )m s 



(38) 



where Ff^ 7T (q 2 ) is B — n transition form factor which is estimated at the q 2 scale, fx, and 
f n are decay constants of K and n meson, respectively. X v i, Y v / represent the rate of 
the uu + dd and ss component in the rj' meson, respectively. We use F^ 7 * (m 2 K ) = 0.35, 
fx = 0.16 GeV, /tt = 0.13 GeV, X v / = 0.57 and Y v > = 0.82. Since rj' contains a small 
uu component, {rj'K°\Qi(2)\B°) has a non zero value. Thus, the B — > rj'Ks process 
has a small contribution from a color suppressed tree diagram. Matrix elements for the 
operator Qi are given by 



(39) 



where the minus sign is a consequence of the parity difference of the initial and final 
states. 

By using the above formulae, numerical values of the matrix elements are found as 
follows: 





Qi 


Q 2 


Qs 


Qa 


Q 5 


Qs 




B -> <pK s 








2.89 


2.89 


2.17 


0.72 


-0.068 


B - V 'K S 


-0.506 


-0.169 


-3.17 


-2.90 


1.55 


-1.96 


0.135 



B Wilson coefficients 

In this Appendix, we give the expressions of the high energy (/j^ ~ 0(mw)) Wilson 
coefficients and of the evolution matrices that are used to calculate the low energy (//& ~ 
0(mb)) Wilson coefficients from the high energy ones. 

Weak scale Wilson coefficients Ci(fiw) are constructed from the SM contributions 
Cf M (fj,w), charged Higgs contributions C^fj-w), gluino contributions Cf(fXw), and chargino 
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contributions Cf(fiw) as 

s-i _ r<SM r< — nSM 

r< — _i_ — /^sm i ^rg 

^3 — ^3 "T ( - / 4 — ( - / 4 i ^4, 

_ /-irSM _|_ sig /~i _ ^iSM _|_ ^ 

C77 = C77 + C77 + C77 + C77' ^89 = Cgg + ^85 + ^89 + ^8g- (40) 

Here, we consider only the SM contributions for the non FCNC process C\{2) ■ Moreover, 
the chargino contributions are taken into account for only CV 7 and Cs g in which the SM 
contributions are suppressed by chirality suppression. O(fiw) scale Wilson coefficients 
are given as follows. [H] 1121 S3] 

Standard Model contributions: 

c sm = 14a. c sm = x _ no. s M = 

1 16vr ' 2 24tt ' 3 24vr v ; ' 

C7f w = C^ 1 = -^E(x t ), Cl M = ^E(x t ), 

C™ = -XtiFfa) + §F 2 (s t )), C£ g M = (41) 

where x t = mf/m 2 ^, 

„, s x(x 2 + llx - 18) -9x 2 + 16s-4, 
E ( x ) = T7T7 7^T^ + T, 777 logo:, 



FAx) 



12(x-l) 3 6(x 
x 3 — 6x 2 + 3x + 2 + 6x log x 



12(x - 1 



, . 2x 3 + 3x 2 — 6x + 1 — 6x 2 logx , . 

2[X) = 12fx-l) 4 • ( > 



Charged Higgs contributions: 



C ?y = -?(&N + f 2N) cot 2 /?+ \f,{x h ) + F 4 (x h )), 



J 7~f — 2\\3 1KHJ ^ 2KHJ ) 3 
2 

where xh = ^tl^Hi 171 h is the charged Higgs mass 



C Z = cot 2 /3 + F 3 (x H )), (43) 



^ x 2 — 4x + 3 + 2x log x , N x 2 — 1 — 2x log x , . 

F3(X) = 27x^TP ' F4(X) = 2(x-l)3 • (44) 
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Gluino contributions: 



ci 



Cf 



r<9 
o 77 



-v/2 



a: 



4G F A t M| 



(E 



p*kspkb -p*ls pis 
1 DL 1 DL 1 DL DL 



1 5 



Y. T D S A[-^C 1 {x kg ) + l -C 2 {x kg )]), 



-V2 



at 



AG F X t Ml 



E r 

kl 



*kspkb p*ls pis 
DL DL DL DL 



7 1 

-^B^Xkg^lg) + -B 2 ( y X k g,X l g) 



E T D S L r DL [lCl(x kg ) - 3 -C 2 (x kg )] ) , 



-V2 



4G F A t Mo 



E r 

kl 



*kspkb p*ls pis 
DL DL DR DR 



^B^X^Xig) + ^B^Xkg.Xtg)^ 



E^r^-^C^) + \C 2 [X k g) 



-V2 



at 



4G F X f M| 



E r 



*fcspfcb p*Zs pis 
DL L DL DR DR 



1 \ r p*fcspfc6 
+ 2_^l nr.l 



fe,; 

) 

DL- 1 - DL 



2 7 1 



[^Ci(x fe9 ) - ^C 2 (a; fe ,%) ), 



-v/2 



a. s 7r 



4G F A t M| 



DL 1 DL 



1 DL 1 DR 



^(-5) E ^(KF 2M - ^«f iM ), 

-V2 a s 



dk 



Ms 
m b 



■g D 2(x k g) 



AG F X t m: 



fE( 



p*kspkb 
1 DL 1 DL 



1 3 

--D^Xkg) + -^3(^9) 



a s n 



V2 



E 



2G jP A t " m' 



pkb p*ks 
1 DL 1 DL 



3 1 
~2 Fl ( X 9k) ~ Q F 2(Xgk) 



-^DR^Dl\-\F,{x gk ) - l -F A { Xgk ) 



(45) 



where x kg =x~ gk 1 = m| /M|, M 3 is the gluino mass, m| represents a mass eigen value of 
the down-type squarks. We define a unitary matrix of the down-type squark mass matrix 

r D = (r D L, r Di? ) as 



(r I) m'-r t 



JJ- D )kl 



m\5 kh 



d 



L(R) 



4 

DL(R) 



dk, 



(46) 



where r is 6 x 3 matrix and m| is the scalar down mass matrix at the weak scale, 
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which can be written as 

_ ., !m 2 ~ L - (i - | sin 2 e w )m 2 z cos 2(3 + m\ -m d (A* d + fj, tan (3) . _ 

'"' / ~~ ' -m d (A d + /i* tan /3) ra| - | sin 2 w m| cos 2/3 + m\ ' ' ^ ' 

The functions appeared above are given by 

-\z 2 dz 



Bi(x,y) 



Ci(x) 
C 2 (x) 

£> 2 (z) 

L> 3 (z) 



o (* + l) 2 ^ + x)(z + y) 

x 2 log a; y 2 log y 



B 2 (x,y) = 



4 ( x - y )(x - 1)2 4(y - x)(y - l) 2 4(x - l)(y - 1) ' 



o (z + l) 2 (z + x)(z + y) 

a; log x |/ log |/ 1 



(x -y)(x - l) 2 (y - x)(y - l) 2 (x-l)(y-l)' 
2x 3 - 9x 2 + 18x - 11 - 6 log x 



36(1 - 


-x) 4 


-16x 3 + 45x 2 - 


36x + 7 + 6x 2 (2x-3) log a; 




36(1 -x) 4 


-x 3 + 6x 2 - 3x 


— 2 — 6a; log x 


6(1 - 


x y 


-x 2 + 1 + 2x lot 


IX 


(x- l) 3 


1 


2x 3 + 3x 2 — 6x - 


f 1 — 6a; 2 log x 



6(1 -x) 4 



„ , . — 3x 2 + 4x — 1 + 2x 2 log x . 

= <— ^ -■ («) 



Chargino contributions: 

j k " 



= ^p^EE^((^"^ 

F 1 j k «fe 

m~- 

-HStiCTj? - H^-^F^J), (49) 

•"b J 

where x--~ = m 2 _/m| fc , and aw is defined from SU{2) L gauge coupling, g, as aw = 

j Xj 

g 2 /(An). m~- and m i( . are mass eigen values of charginos and scalar up quarks, respec- 
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tively. Guli Hul, and Hjjr are defined as 

r~ijki t r* -pki 

^UL ~ v jl L ULi 

HtR = V* 2 (T UR Y U V C 



H$t = U j2 (T UL Y d r/g 



ki 



CKM 



\ ki 



/g, 



(50) 



where IV = (Tm, IV^) is a unitary matrix of the scalar up mass matrix which is defined 
in the same manner as the scalar down mass matrix, and V and U are diagonalizing 
matrices of the chargino mass matrix which are defined as 



(UmcV^ij = m^Sij, m c 



M 2 \/2mw sin ft 

\/2mw cos ft fi 



(51) 



Next, we have to estimate the Wilson coefficients at the scale O(fib) from the Wilson 
coefficients at the scale 0(fj,w) using the renormalization group equations [21]. To this 
end, we can use the evolution matrix as follows. 



1,...6) 



(52) 



Here, we take account of only the leading order QCD corrections for simplicity. At the 
leading order, the evolution matrix U(fib, ^w) is estimated as 



w j 



- Oi s (/jLb) 



JO) 

23() 



D 



(53) 



where ft = -jNc — \f is the coefficient of the one loop ft function of the QCD gauge 
coupling. Here, Nc and / are number of color and number of active flavor at the scale /i 
(/it, < /i < Hw)-, respectively. 7^ = V~ 1 ^°"> T V is the diagonalized matrix of 7^ which 
is given by 

6 \ 

2-22 



7 



(0) 



Using N c = 3,f 



U{/ib,Hw) 



/-2 
6 




\o 

5, a(fib) : 

/ 1.117 
-0.267 
-0.0014 
0.0034 
-0.0010 

\ 0.0039 



-2 — 

9 

-22 










6 _ U 
u 9 



-2/ 
9 



I 

-2 + ¥ 



2/ 
3 





-2 

-% 
9 

2 

zM. 

9 



t 

3 

-6 
-16 + 



2/ 



(54) 



0.22, a{fi w ) = 0.1176 [Sj we obtain 



-0.267 
1.117 

0.0120 
-0.0273 

0.0079 
-0.0341 






1.134 
-0.305 
0.0107 
-0.0484 






-0.210 
0.987 
0.0376 
-0.162 






0.0079 






0.0791 



-0.0197 0.0541 
0.928 0.0541 



For C7 7 (mb) and C% g {mb) 



0.309 

we can use the following formulae [2j^ 



(55) 



1.692 / 



C 7l {m b ) = 0.695CV 7 (/in/) + 0.085C 8a (// w ) 
C 8g {m b ) = 0.727C 8g {»w) - 0.074C 2 (/i W ). 



0.158C 2 (/x 



w 1 



(56) 
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